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Abstract 

>: 

l/*~) , We consider the radially vibrating spherical quantum billiard as a represen- 

1/"") ' tative example of vibrating quantum billiards. We derive necessary conditions 

for quantum chaos in d-term superposition states. These conditions are sym- 
metry relations corresponding to the relative quantum numbers of eigenstates 
considered pairwise. In this discussion, we give special attention to eigenstates 
f~^ ■ with null angular momentum (for which the aforementioned conditions are au- 

""^ \ tomatically satisfied). When these necessary conditions are met, we observe 

numerically that there always exist parameter values for which the billiard 
behaves chaotically. We focus our numerical studies on the ground and first 
J* ■ excited states of the radially vibrating spherical quantum billiard with null an- 

. ,_h [ gular momentum eigenstates. We observe chaotic behavior in this configuration 

^^ ■ and thereby dispel the common belief that one must pass to the semiclassical 

J-j \ (h — ► 0) or high quantum number limits in order to meaningfully discuss quan- 

tum chaos. The results in the present paper are also of practical import, as the 
radially vibrating spherical quantum billiard may be used as a model for the 
quantum dot nanostructure, the Fermi accelerating sphere, and intra-nuclear 
particle behavior. 

1 Introduction 

There has been considerable research in the last twenty years that seeks to marry 
quantum mechanics and dynamical systems theory into a coherent whole, fl 0] In 
particular, the concept of quantum chaos extends the notions of classical Hamil- 



o 



toman chaos to the quantum regime. There are three types of quantum chaotic 
behavior: "quantized chaos" ( "quantum chaology" ) , "semiquantum chaos," and 
genuine "quantum chaos." Quantum chaology concerns the quantum structure 
of classically chaotic systems, semiquantum chaos refers to the chaotic dynamics 
of coupled classical and quantum systems, and genuine quantum chaos refers 
to chaotic dynamics of fully quantum systems. || No example of the third type 
of quantum chaos has been established, so the existence of such systems is an 
open question. 

In the present paper, we discuss semiquantum chaos in the context of vi- 
brating billiard systems. We review the recent results of Liboff and Porter H 
and discuss them in further detail. We treat a two-term Galerkin projection 
(superposition state) of the radially vibrating sphere and prove that only 2- 
term superpositions whose normal modes have common rotational symmetry 
behave chaotically. We extend this theorem to arbitrary superposition states 
by applying it pairwise. In the proof of this theorem, we establish integrable 
(non-chaotic) behavior by showing that the evolution equations reduce to a 
two-dimensional autonomous dynamical system, whose non-chaotic properties 
are known. [Q We then discuss two examples: one integrable superposition and 
one chaotic one. We compute Poincare maps and thereby reveal chaotic char- 
acteristics such as regions of ergodicity and KAM islands. 

We also discuss the present results with respect to the phenomenology of 
quantum chaos. The chaotic behavior in the radius a and conjugate momentum 
P corresponds to classical Hamiltonian chaos. The normal modes ip n im depend 
on the radius, so they exhibit quantum-mechanical wave chaos. We also observe 
chaos in the Bloch variables (x, y, z), which correspond to quantum-mechanical 
probabilities. The dynamical equations describing the present system corre- 
spond to a two degree-of- freedom (dof) Hamiltonian system, where one degree- 
of-freedom is classical (corresponding to the so-called one "degree-of- vibration" 
(dov) radial motion) and one is quantum- mechanical (corresponding to the cou- 
pling coefficient fx). By coupling a single classical dof (which is necessarily 
integrable) with a single quantum dof (which must also be integrable), we ob- 
tain a genuinely chaotic system that provides an example of semiquantum chaos 
since it consists of a classical system coupled to a quantum one. We remark that 
we do not need to pass to the semiclassical (H — ► 0) or high quantum num- 
ber limits in order to observe chaos, as is commonly considered requisite for a 
meaningful analysis of quantum chaos. H 

The radially vibrating spherical quantum billiard has several practical appli- 
cations that complement its theoretical import. The most important one is that 



it may be used as a model for the quantum dot nanostructure.|10 At low tem- 
peratures, this microdevice component experiences vibrations due to zero-point 
motions, and at higher temperatures, it exhibits vibrations due to natural fluc- 
tuations. Another application is that the radially vibrating spherical quantum 
billiard generalizes Fermi's "bouncing-ball model" of cosmic ray acceleration, pj 
Additionally, the radially vibrating spherical quantum billiard models the in- 
tradynamics of the nucleus, as the 'liquid drop' and 'collective' models of the 
nucleus include boundary vibrations. Consequently, the importance of the ra- 



dially vibrating spherical quantum billiard lies not only in its expansion of the 
theory of quantum chaos but also in its applicability to problems in nuclear, 
atomic, and mesoscopic physics. 

2 Statement of the Problem 

The spherical quantum billiard addresses the quantum dynamics of a particle 
of mass too confined to the interior of a spherical cavity of mass M ^S> too with 
smooth walls of radius a. The radius vibrates in an a priori unspecified manner, 
so that a = a(t). A two-component superposition state (Galerkin projection) of 
this quantum billiard is given using Dirac notation by 

|V>(r,0,<£,t;a(t))) = Ai(t)\nlm,t) + A 2 (t)\ril'm',t), (1) 

where A\{t) and ^(i) are complex amplitudes. The numbers {n, l,m} are, re- 
spectively, the principal, orbital, and azimuthal quantum numbers. The eigen- 
states of the present system are products of spherical Bessel functions and spher- 
ical harmonics. H In coordinate representation, 



(r\nlm,t) =ip n im(r,6,(f>,t;a(t)) = * I —rrr; I -, r ) ji I — rr ) Y lm (d,(f>), 

V am \3i+i{xi n )J \a{t) ) 

(2) 

where xi n is the nth zero of ji, the spherical Bessel function of order I. 

For the system at hand, the time-dependent Schrodinger equation is given 

by 

ih^ =Kil> ) r< a(t), (3) 

where the quantum-mechanical Hamiltonian K, the kinetic energy of the parti- 
cle, is 

fi 2 
K = --^V 2 . (4) 

2to V ' 

The total Hamiltonian of the system is 

where P is the momentum of the billiard boundary, and V = V(a) is the 
potential of the billiard surface. The potential energy V and kinetic energy 
P 2 /2M of the billiard walls are classical quantities, and the confined particle 
is quantum-mechanical. For this semiquantum system, we utilize the Born- 
Oppenheimer approximation [|2|, so that only the particle kinetic energy K is 
inserted into the Schrodinger equation (pf). In this adiabatic approximation, 



which is commonly used in mesoscopic physics, we are ignoring the effects of 
Berry phase. [fl3[ 

Taking the expectation of (pf) using the superposition state (pi) gives 
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V V) = — [eil^il 2 + e 2 |A 2 | 2 ] = K(Ai,Aa,a), 



iH{ip 



dip 
~dt 



ifr 



A X A\ + A 2 A* 2 + ^nl^il 2 + v 12 A x Al + v 21 A 2 Al + v 22 \A 2 \ 



where the energies of the two terms are given by 

fe2 2 
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(7) 



3 Integrable Configuration 



Examining the superposition of 1 100) and 1 110) using (pi) and orthogonality of 
spherical harmonics shows that v\\ = v\ 2 = v 2 \ = v 22 = 0. (Note that V\\ 
and v 22 vanish no matter which eigenstates one considers.) Equating the inner 
products ffl) of both sides of the Schrodinger equation gives equations of motion 
for the complex amplitudes: 



%A\ = r-^eiAi, 



iA 2 



h^ e2M > 



(8) 



which are integrated to yield 



At (t) = C\ exp 



a~ 2 (t)dt 



, A 2 {t) = C 2 exp 



-^ / a~ 2 (t)dt 



(9) 



From (ph, one obtains a Hamiltonian in the radius a and conjugate momentum 
P: 



H 



2M 



+ K(A 1 ,A 2 ,a) + V(a) 



2M 



[e.lC.l 2 + e 2 \C 2 \ 2 ] 



V{a). (10) 



A Hamiltonian with no explicit time-dependence and one dof corresponds to a 
two-dimensional autonomous system, which is known to be non-chaotic. pi |l2[ 
When there are no coupling terms, the degree-of-freedom of the resulting Hamil- 
tonian corresponds to the degree-of- vibration of the quantum billiard, which is a 
measure of the number of distance dimensions that undergo oscillations. When 
a two-term superposition has a non-vanishing coupling coefficient, the number 
of degrees-of-freedom of the resulting Hamiltonian system is equal to the num- 
ber of dov of the billiard plus one. In particular, this means that a superposition 
state of a quantum billiard with more than one dov (such as the two dov rect- 
angular quantum billiard) is expected to behave chaotically even if every one of 
its coupling coefficients vanishes. 



Hamilton's equations for the present integrable configuration are 

P dH ■ dV A dH 

M oP aa a A da 

where the energy parameter A is given by 

A = 2(e 1 |Ci| 2 + e 1 |C 2 | 2 ) > 0. (12) 

The bifurcation structure of ( fll| ) has been studied for quartic potentials V(a).jfl| 

4 Necessary Conditions for Chaos in k Coupled 
States 

Consider the superposition 

tp = Aii> qi + A 2 ijj q2 + ■■■ + A k ij) qk , (13) 

where qi = (m, U, rrii) is a vector of quantum numbers. If there does not exist 
a pair of normal modes in the fc-state superposition (U3) with common angular 
momentum quantum numbers (i.e., there is no pair {i,i'} such that h = h> 
and rrii — my ) , then inserting (O) into the Schrodinger equation (0) returns a 
diagonal quadratic form 

MA\ + ■■■ + A k A* k = vn\A x \ 2 + ■■■ + v kk \A k \ 2 , (14) 

as all the cross terms UijAiA* have vanishing coupling coefficients i/y by orthog- 
onality of spherical harmonics with different angular momenta. The diagonal 
terms in (H4) stem from the Laplacian. As above, we obtain the Hamiltonian 



H(a,P) = ^ + ^Y.£ t \a\ 2 + V(a), (15) 



i=i 



where the Ci are constants. The superposition ( |13[ ) is non-chaotic, because the 
Hamiltonian ( |l5| ) is autonomous with one dof. 

We thus conclude that a necessary condition for chaotic behavior in an ar- 
bitrary finite superposition state of the radially vibrating spherical quantum 
billiard is that at least one pair of normal modes in the eigenfunction expan- 
sion have common angular momentum quantum numbers. In particular, by 
considering small n, and 7ty , we obtain a chaotic superposition for eigenstates 
with small energies. This even holds for some superpositions that include the 
ground state! In most studies of quantum chaos, one must take the semiclas- 
sical {h — ► 0) or high quantum- number limits in order to meaningfully study 
quantum chaos. 0, pi In such studies, termed "quantum chaology," one considers 
the quantum signatures of classically chaotic systems in these regimes. In the 
present system, on the other hand, we obtain genuinely chaotic behavior in a 
semiquantum system. This phenomenon is often called "semiquantum chaos." B 



5 Chaotic Configuration 

As an example of a chaotic configuration of the radially vibrating spherical 
quantum billiard, consider the superposition state consisting of the ground and 
first excited states with null angular momentum 

\ip(n,l,m))=A 1 \100) + A 2 \200) ) (16) 

which gives the wavefunction 



4>(r,t) = A^mMr,^ 1 ^ 1 + A 2 {t)a 2 i> 2 {r,t)e- lE ^ , (17) 



where 



, . / rurr\ . . sin(x) V2 

Wn(r,t)= jo \-7-r ) , Jo s = , a n = -g— - — -. (18) 

\a{t)J x a*ji(rm) 

The superposition (19) has a coupling coefficient p = p\ 2 = 4/3. 

Equating coefficients in the quadratic form (p) gives the matrix equation 

2 

iA n = Y,D n kA k , (19) 

fe=l 

where D = (Dij) is the Hermitian matrix 

ei 

ha ■ 

and the energy coefficient e-j is given by 

_ CM) 2 



D=[Ea T a ) > (20) 



2mo 



ie{i,2}. (21) 



Defining the density matrix[g| by p 9 „ = AqA^, introducing (dimensionless) 
Bloch variables x = p\ 2 + p 2 \, y = i{p 2 \ — P12), and z = p 22 — pu, and using 
(|19|), we obtain the following three differential equations: 

. _ uj y 2pPz . _ uj x . _ 2pPx 

ar Ma a 1 Ma 

In these equations, u>q = (e 2 — e\)/h. Rewriting the kinetic energy K{A\, A 2 , a) 
in terms of the Bloch variable z gives 

K{z 1 a) = ^(e+ + ze-) 1 e ± = -(e 2 ± ei ). (23) 

Inserting K{z, a) into the Hamiltonian (0) gives Hamilton's equations: 

P ■ dV 2 

°=M< P = -lTa + -^ + + e - {Z -» X)] - (24) 



Equations (122) and ( |24| ) constitute a set of five coupled nonlinear ordinary 
differential equations, which can be shown to be equivalent to a two degree-of- 
freedom Hamiltonian system. The constants of motion of the present system 
are the radius of the Bloch sphere 

x 2 +y 2 + z 2 = \A 1 \ 2 + \A 2 \ 2 = l (25) 

and the energy (total Hamiltonian) 

H=^- [ +V(a)+K{z,a). (26) 

The equilibria of equations ( p2| . p3) satisfy x = y = 0, z = ±1, P = 0, and 
a — a±, where a± satisfies the equation 

oa a J 

where the subscript of a± corresponds to the sign of z = ±1. Assuming that 
V(a) + K(z,a) has a single minimum in a, these equilibria are elliptic. |)|, |ll| 
(That is, every eigenvalue of the Jacobian of the linearized system is purely 
imaginary.) For the harmonic potential 

V(a) = ^(a-a ) 2 , (28) 

the total energy of the billiard's boundary is given by 

V(a) + K(a) = \a - a ) 2 + '+ - ' Z€ ~ . (29) 

With this choice of potential, equation (p^) becomes 

«-„„ = f± ,30) 

The solutions of ([30|) for the equilibrium radii a± correspond to the e± values. 
One computes that a+ > a- > aoj from which it follows that a_ < a(0) < a + , 
so a(t) remains bounded in the interval [a_, a+].B 

Now consider the superposition of the first k null angular momentum eigen- 
states, 

k 
^W(r-) = ^ A n (t)a n Mr,t). (31) 

71=1 

In order to analyze this configuration, one first examines the 2-term superposi- 
tion 

'ipnq = A n (t)Onipn(r, t) + A q (t)a q ip q (r, t), n < q, (32) 



and then superposes the couplings one obtains from each ip nq as n and q run 
from 1 to k in order to obtain dynamical equations for the amplitudes A;. One 
computes the coupling coefficients fi ng to be 

Mn<2 = 2— 7-r- ■ r, n ^ q. (33) 

a{t)(n + q)(q - n) 

The dynamical equations for A+ are described by a k x k matrix and are a 
straightforward generalization of ( |l9| , f20|) . 

5.1 Numerical simulations 

The analysis for 2-term superpositions of null angular momentum eigenstates 
follows that for the general case.|9j In the present case, the necessary conditions 
for chaotic behavior are satisfied automatically, because the quantum numbers 
m and I vanish for every normal mode under consideration. Consequently, any 
fc-term superposition (k > 2) of null angular momentum eigenstates exhibits 
chaotic behavior. We consider numerical simulations for the coupling of the 
ground state and first excited state of a billiard residing in a harmonic potential. 
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Figure 1: Poincare Section (x = 0) in the (a,P)-plane illustrating that not all 
invariant tori are destroyed in the present configruation. 

Figure |l| shows a Poincare map in the (a, P)-plane corresponding to x = 0, 
and Figure^ snows a Poincare section projected onto the (x, y)-plane for P = 0. 
For each of these two plots, we used the parameter values U = 1, M = 10, 
m = 1, Vo/oq = 5, and clq = 1.25. The initial conditions for the two figures 
are x(0) = sin(0.957r) « 0.156434, y(0) = 0, z(0) = cos(0.95tt) w -0.987688, 
a(0) w 1.6, and P(0) w 9.45. 

The chaotic behavior of this configuration is evident in both plots, although 
there is clearly still some non-chaotic structure present. In the language of 
KAM theory, some of the nonresonant tori persist for the present choice of 
initial conditions fa, 121. One may also choose initial conditions corresponding 
to a different level of persistence of the resonant tori. For example, Figure pi 
shows an x = Poincare map in the (a, P)-plane with the same initial conditions 
and parameter values as above, except a(0) = 3 and P(0) = 10. Figure shows 
a P = Poincare map in the (x, y)-plane for these conditions. There are fewer 
invariant tori in these two figures than there are in Figures wfl- 









Figure 2: Poincare Section (P = 0) of the Bloch sphere projected onto the 
(x, y)-plane. The structure in this diagram likewise illustrates the survival of 
some invariant tori. 

6 Phenomenology 

In contrast to the present quantum-mechanical context, we note that for the clas- 
sical radially vibrating spherical billiard, every orbit with null angular momen- 
tum is well-defined and invariant under radial oscillations of the boundary. H 
Due to conservation of angular momentum, one finds that for the stationary 
spherical classical billiard, the enclosed particle sweeps out an annular domain 
of constant inner radius. || Vibration of the wall of the sphere destroys this 
constant, and chaotic motion is expected to develop. In the present quantum- 
mechanical context, we note that null angular momentum wavefunctions are 
composed only of spherical waves. The nodal surfaces of these wavefunctions 
are likewise spherical. Accordingly, the chaotic signature of this configuration 
in real space is the sequence of intersections with a fixed radius that nodal sur- 
faces make at any instant subsequent to a number of transversal times. H This 
latter condition is consistent with the standard long-time behavior of chaotic 
dynamical systems. || 

In the language of Blumel and Reinhardt fl , vibrating quantum billiards are 
an example of semiquantum chaos. One has a classical system (the walls of the 
billiard) coupled to a quantum-mechanical one (the enclosed particle) . Consid- 
ered individually, each of these subsystems is integrable, as each contributes a 
single dof. When they are coupled, however, one observes chaotic behavior in 
both of them. The classical variables (a, P) exhibit Hamiltonian chaos, whereas 
the quantum subsystem [x, y, z) is truly quantum chaotic. Chaos on the Bloch 
sphere is an example of quantum chaos because the Bloch variables (x, y, z) cor- 
respond to the quantum probabilities of the wavefunction. Additionally, each 
individual normal mode ip n depends on the radius a(t), so each eigenfunction is 
an example of quantum-mechanical wave chaos for chaotic configurations of the 
billiard. Moreover, because the evolution of the probabilities \Ai\ 2 is chaotic, the 
wavefunction ifj in the present configuration is a chaotic combination of chaotic 
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Figure 3: Poincare Section (x = 0) in the (a, P)-plane for slightly different 
initial conditions in which fewer invariant tori persist. This is in accord with 
KAM theory. 



normal modes. Finally, we note that if we quantize the motion of the billiard 
walls, we would obtain a higher-dimensional, fully-quantized system that ex- 
hibits so-called quantized chaos. B In particular, the fully quantized version of 
the present system would require passage to the semiclassical limit in order to 
observe quantum signatures of classical chaos. 

7 Conclusions 

We considered the radially vibrating spherical quantum billiard as a represen- 
tative example of vibrating quantum billiards. We derived necessary conditions 
for quantum chaos in /c-term superposition states. We gave special attention to 
eigenstates with null angular momentum, for which these conditions are auto- 
matically satisfied. We examined a numerical simulation of the superposition of 
the ground and first null angular momentum excited states. We observed chaotic 
behavior in this configuration, thereby dispelling the common belief that one is 
required to pass to the semiclassical {% — ► 0) or high quantum number limits 
in order to meaningfully study quantum chaos. 
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Figure 4: Poincare Section (P = 0) of the Bloch sphere projected onto the 
(x, y)-plane. The initial conditions in this plot are the same as those in Figure 
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